Roughening and inclination of competition interfaces 
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We study the competition interface between two clusters growing over a random vacant sector of 
the plane in a simple set-up which allows us to perform formal computations and obtain analitical 
solutions. We demonstrate that a phase transition occurs for the asymptotic inclination of this 
interface when the final macroscopic shape goes from curved to non curved. In the first case it 
is random while in the second one it is deterministic. We also show that the flat case (stationary 
growth) is a critical point for the fluctuations: for curved and flat flnal proflles the fluctuations 
are in the KPZ scale (2/3); for non curve final profile the fluctuations are in the same scale of the 
fluctuations of the initial conditions, which in our model are Gaussian (1/2). 

PACS numbers: 64.60.Ak 64.60.Ht 
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Introduction The behavior of the interface of a grow- 
ing material has been investigated using the Eden model 
ballistic deposition and other random systems. Typ- 
ically, the growing region converges to an asymptotic de- 
terministic shape and its fluctuations depend on the ge- 
ometry of the initial condition ■ ^ l^ss well stud- 
ied phenomenon is the competing growth of two mate- 
rials. The interface between two growing clusters (com- 
petition interface) presents a random direction on the 
same scale as the deterministic shape 0, ^2 ^3 ■ In this 
letter we describe quite explicitly this phenomenon in a 
simple model. On gronds of universality, this will pro- 
vide a guide to understand the interplay between the 
asymptotics of the competition interface and the final 
macroscopic shape in models with different growth and 
competition mechanisms. 

We determine the inclination of the competition interface 
for a growth model called "last passage percolation" in a 
random sector of the plane of angle 0. The growth inter- 
faces are mapped into particle configurations of the to- 
tally asymmetric simple exclusion process in one dimen- 
sion (TASEP) [3. Under Euler space-time rescaling, the 
particle density of the TASEP converges to a solution of 
the Burgers equation. This equation has travelling wave 
solutions (shocks) corresponding to the case 9 > 180°, 
and rarefaction fronts corresponding to 9 < 180°. A per- 
turbation at one site of the initial particle configuration 
(called a second class particle) follows a characteristic of 
the equation or the path of a shock. To establish our 
results, we map the competition interface linearly onto 
the path of the second class particle. 

The growth model The random sector is 
parametrized by the asymptotic slope of its sides. 
Let A G (0, 1] and p G [0, 1) and define a random path 
70 = ho{j))jez C Z2 with 71 (0) ^ (1,0), 7o(0) - (1, 1), 



7-i(0) = (0,1) as follows. Starting from (0,1), walk 
one unit up with probability A and one unit left with 
probability 1 — A, repeatedly, to obtain 7g = (7o(j))i<o- 
Then, starting from (1,0) walk down with probability p 
and right with probability 1 — p to get 7q = (7o(j))i>o- 
7q has asymptotic orientation (A — 1,A) while 7q has 
asymptotic orientation (1 — p, —p). Let Co be the 
sector with boundary 70, containing the first quadrant; 
its asymptotic angle 9 = 9x,p e [90°, 270°) is the 
angle between (A — 1,A) and (1 — p,—p). Notice that 
9 e [90°, 180°) if and only if p < A. 

The path 70 is the growth interface at time 0. The dy- 
namics are then defined as follows. For each z G Co and 
each t > 0, we have a label (Tt{z) G {0,1,2}. The la- 
bel is if z is unoccupied at time t, and is 1 or 2 if z 
belongs to cluster 1 or 2 respectively. Once occupied, a 
site remains occupied and keeps the same value forever. 
Initially, set ao{z) = 1 for all z G 7o, o-o(z) = 2 for all 
2: G 7o and ao{z) = for all z G Co\7o- Independently 
each vacant site z G Co\7o becomes occupied with rate 
1 provided z — (1,0) and z — (0,1) are occupied. Let 
G{z) be the time at which site z becomes occupied. At 
this time crt(z) assumes the value crt(z) where z is the 
argument that maximizes G(z — (1, 0)) and G(z — (0, 1)). 
Thus when a site becomes occupied it joins the cluster of 
whichever of its two neighbours (below and to the left) 
became occupied more recently. The label of the site 

(1.1) may be left ambiguous, but we stipulate that site 

(1.2) always joins cluster 1, and site (2, 1) always joins 
cluster 2. 

The process {Gj, G^), where is the set of sites z G Co 
such that (Tt(z) = fc, describes the competing spatial 
growth model. The growth interface at time t is the 
polygonal path 74 composed of sites z G Co such that 
G(z) < t and G(z -f (1,1)) > t. The competition inter- 
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FIG. 1: Growth and competition interfaces. 



face (fi = ((y3„)N is defined by (po = (I7 1) and, for n > 0, 
V'n+i = (fin + (1,0) if (fn + (1,1) £ G;^ and ipn+l = 
ipn + (0, 1) if (fin + (1, 1) S G^. Note that ip chooses lo- 
cally the shorter step to go up or right, so that ipn+i is the 
argument that minimizes {G{ipn + {1, 0)),G{(pn + {0, 1))}. 
This competition interface represents the boundary be- 
tween those sites which join cluster 1 and those joining 
cluster 2 (see FigureP). The process V'(i) = {I{t),J{t)) 
defined by i(j{t) = ipn for t S [G{ipn),G{(pn+i)) gives the 
position of the last intersecting point between the com- 
petition interface (p and the growth interface 74. 

In we prove that with probability one. 
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where a £ [0, 90°] is given by 
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and J7 is a random variable uniformly distributed in [1 — 
2A, l-2p]. 

Simple exclusion and second-class particles The 

totally asymmetric simple exclusion process {rjt, t > 0) 
is a Markov process in the state space {0, 1}^ whose el- 
ements are particle configurations. ??t(j) = 1 indicates a 
particle at site j at time t, otherwise ijt{j) = (a hole 
is at site j at time t). With rate 1, if there is a particle 
at site j, it attempts to jump to site j + 1] if there is a 
hole at j + 1 the jump occurs, otherwise nothing happens. 
The basic coupling between two exclusion processes with 
initial configurations 770 and rjQ is the joint realization 
{rit,ri[) obtained by using the same potential jump times 
at each site for the two different initial conditions. Let 770 
and rjQ be configurations of particles differing only at site 
X(0) = 0. With the basic coupling, the configurations at 
time t differ only at a single site X{t), the position of a 
so-called second-class particle. Such a particle jumps one 
step to its right to an empty site with rate 1, and jumps 
backwards one step with rate 1 when a (first class) par- 
ticle jumps over it. 



If r;o is distributed according to the Bernoulli product 
measure with density A for j < and p for j > Q, then 
the asymptotic behavior of X{t) shows a phase transition 
in the line X = p: with probability one. 
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where [/ is a random variable uniformly distributed in 
[1 — 2A, 1 — 2p\ (0,0,^3 fo'^ deterministic case and 
0,11, 13, El for the random case). 

The limits (jSJ are based on the following hydrodynamic 
limits. If 770 is distributed with the product measure with 
densities A and p as before, then the macroscopic density 
evolution is governed by the Burgers equation: 

^^'^^^Y.fi^^)'lt/e{x)= / f{r)u{r,t)dr (4) 

with probability one for all / : R ^ K with compact sup- 
port, where u(r, t) is the solution of the Burgers equation 

+ - ^))) = 0' reR,t>0 

with initial condition u{r,0) = A for ?' < and p for 
r > 0. If A = p the solution is constant, if A < p it is a 
shock: 



u{r,t) 



A if r < (1 - A - p)i 
p if r > (1 - A - p)i 



(5) 



and it is a rarefaction front if A > p: 



u(r, t) 
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if r < (1 - 2\)t 
2(X^ if(l-2A)t<r<(l 
if r > (1 - 2p)t 



2p)t 
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(^3,01 for initial product measures and p^ J I ^ for ini 



tial measures satisfying with t = 0; also [lll2l) 



The characteristics v{a,t), corresponding to the Burg- 
ers equation and emanating from a, are the solutions of 
dv/dt = 1 — 2u{v,t) with v{Q) = a. The solutions are 
constant along the characteristics. When two character- 
istics carrying a different solution meet, they give rise to 
a shock. There is only one characteristic emanating from 
locations where the initial data is locally constant and 
there are infinitely many characteristics when there is a 
decreasing discontinuity. In particular, if the initial con- 
dition is u{r, 0) = A for r < and u(r, 0) = p for r > 0, 
then the characteristics Vr{t) emanating from the point 
Vr{0) = r are given by Vr{t) = r + (1 — 2X)t if r < 
and Vr{t) = r -I- (1 — 2p)t if r > 0. For r = there are 
two cases. When A < p, the characteristics emanating 
from positive sites are slower than those emanating from 
negative sites. They collide, giving rise to a shock Q 
traveling at speed 1 — A — p. When A > p there are in- 
finitely many characteristics emanating from the origin: 
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for each s G [1 — 2 A, 1 — 2p\ the hne i'o(i) — st\s a. charac- 
teristic emanating from 0. The hmits lO show that the 
second-class particle follows the characteristic when there 
is only one (that is, when A = p), that it follows the shock 
when the initial condition has an increasing discontinuity 
and that it chooses uniformly one of the characteristics 
emanating from a decreasing discontinuity. 

Growth and simple exclusion Rost relates the 
simple exclusion process to the growth model as follows. 
Consider initial configurations rj^ for the exclusion pro- 
cess in which 770(0) = and ?7o(l) = 1- Elsewhere let 
770 be distributed according to Bernoulli product mea- 
sure with density A for j < and p for j > \. De- 
fine the initial growth interface 70 by 70(0) = (1, 1) 
and 7o(j) - 7o(j -!) = (!- r?o(j), -?7o(j)); then 70 
has the same distribution as before. Label the parti- 
cles sequentially from right to left and the holes from 
left to right, with the convention that the particle at 
site 1 and the hole at site are both labeled 1. Let 
Pj(0) and Hj{0), j S Z be the positions of the parti- 
cles and holes respectively at time 0. The position at 
time t of the jth particle Pj{t) and the ith hole Hi{t) 
are functions of the variables G{z) with z € Co\7o (de- 
fined earlier for the growth model) by the following rule: 
at time G((i,j)), the jth particle and the ith hole in- 
terchange positions. Disregarding labels and defining 
r]t{Pj{t)) = 1, r]t[Hj{t)) = 0, j S Z, the process 77* indeed 
realizes the exclusion dynamics. At time t the particle 
configuration 774 and the growth interface 74 still satisfy 
the same relation as 770 and 70. This connection yields 
the following shape theorem for the growth model. Al- 
most surely, 
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where h{r) — h\^p{r) is related to the hydrodynamic limit 
(TO by h'{r) = w(r,l)/(l -u(r,l)). 

Second class particles and competition interfaces 

A key tool in proving 11I2() is the observation |^ that the 
process given by the difference of the coordinates of the 
competition interface I{t) — J{t) behaves exactly as the 
second class particle initially put at the origin. To see 
this call the particle at site 1 *particle and the hole at 
site *hole, and call this couple *pair. The dynamics 
of the *pair is the following: it jumps to the right when 
the *particle jumps to the right, and it jumps to the left 
when a particle jumps from the left onto the *hole. The 
*pair then behaves as a second class particle. The only 
difference is that it occupies two sites while the second 
class particle occupies only one. The labels of the *par- 
ticle and *hole change with time. At time they both 
have label 1 and the labels of the *pair are represented 
by the point (^0 = (Ijl); the initial value of the com- 
petition interface. If, say, 0(2,1) < G(l,2), then the 
*particle jumps over the second hole before the second 
particle jumps over the *hole (see Figure|21). In this case, 
the labels of the *pair at time G(2, 1) are (2, 1), which is 



exactly the argument that minimizes {G(2, 1), 0(1, 2)}; 
thus, after the first jump of the *pair. its labels are given 
by (pi . By recurrence, (p„ gives exactly the labels of the 
*pair after its nth jump. Therefore the labels of the '''par- 
ticle and *hole are J{t) and /(<), respectively. In addi- 
tion, J{t) — 1 is exactly the number of jumps that the 
*pair has made backwards up to time t, and I{t) — 1 is 
the number of its jumps forwards. This shows that if the 
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FIG. 2: Pair representation of second class particle. 

exclusion and the growth process are realized in the same 
space, X{t) ~ I{t) — J{t). As a consequence of this and 
we get the following behavior for ■0(i) that implies 
(P) and lO. Almost surely. 
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((l-p)(l-A),Ap) ifA<p 
i((C/ + l)2,(t/-l)2) ifA>p 
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where [/ is a random variable uniformly distributed in 
[1-2A, l-2p]. 

To prove JH)) for A > p recall that Pi{t) is the position 
of the 1st particle at time t. Thus J{t) is the number of 
particles that at time zero were to the left of Pi(0) = 1 
and at time t are to the right of X{t). Therefore, J{t) is 
equal to the number of particles between X{t) and Pi{t) 
at time t. By the law of large numbers, Pi{t)/t converges 
to 1 — p and, hy X{t)/t converges to U. Hence J{t)/t 
converges to the integral of the solution of the Burgers 
equation at time 1 (w(r, 1) given by ©) in the interval 
{U, 1 - p). Taking /(r) = l[re[u,i-p]] in 



Jit) 



Pi{t) 
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Vt{j) 



u{r,l)dr^ -{1~U)- 



Analogously, since I{t) is the number of holes to the right 
of Hi{0) = at time zero and to the left of X{t) at time 
t and Hi{t)/t converges to —A almost surely, we obtain 
(O for A > p. For X < p the same argument works by 
substituting U above by 1 — A — p, the limit position of the 
second class particle in this case, and taking the solution 
w(r, 1) given by |(SJ). 

Fluctuations For 9 > 180° the second class particle 
has Gaussian fluctuations produced by the initial pro- 
file 0. This together with the relation above implies 
that under a diffusive scaling {I{t),J{t)) converges to a 
bidimensional Gaussian distribution with a non-diagonal 
covariance matrix computed explicitly 

To understand the fluctuations for 6 < 180° we relate the 
models to a directed polymer model. For each z G Co\7o 
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let = G{z) - max{G(z - (1,0)), G{z - (0, 1))}. Then 
{wz, z £ Co\7o) is a sequence of i.i.d random variables 
with an exponential distribution of mean 1. Let n(z, z') 
be the set of all directed polymers (or up-right paths) 
(zi,...,Zn) connecting z to z' , and let G{z',z) be the 
maximum over all tt G n(z',z) of i(7r), the sum of Wz 
along the polymer tt. Each site z has energy — w^, 
and the polymer tt has energy — i(7r). Thus — G(z',z) 
is the minimal energy, or ground state, between z' and 
z. There exists a unique polymer M(z', z) in Tl{z' , z) 
that attains the maximum. We say that the semi-infinite 
polymer (2:„)n is maximizing if for all ?i < to we have 
(z„,...,Zm) = M(z„,Zm). Every semi-infinite maxi- 
mizing polymer (zn)N has an asymptotic inclination e*" 
0. In the competition model, G(z) = G(7o,z) and for 
k = 1,2, is the set of sites z such that M(7o, z) orig- 
inates from 7g . Denote by the roughening exponent of 
semi-infinite maximizing polymers. 

For Q = 180° (A = p) the process is stationary and 
the connection is explicit. Running the process for- 
ward and backward we extend G(z) to all z G Z^; 
G+ = (G(z), z G 1?) and G" = (-G(-z), z G 1?) are 
indentically distributed. We define the forward competi- 
tion interface starting at z, (p^ = (<^^)n, by setting fa — z 
and putting <p^+i equal to the argument of the minimum 
between G{ipf^ + (1,0)) and G(<p^ -|- (0, 1)), and the back- 
ward semi-infinite polymer starting at z, = (M^)n, 
by setting M§ = z and putting M^^i equal to the ar- 
gument of the maximum between G{Mjl — (1,0)) and 
G(M^ - (0, 1)). Note that (p = (^(^'i' and that Af^ is 
a semi-infinite maximizing polymer. Together with the 
duality relation ip'{G^) ~ M^{G~), this shows that the 
forward competition interface has the same law as the 
backward semi-infinite maximizing polymer and, in par- 
ticular, they have the same fluctuations, so that x = ^■ 

For 6 < 180° (A > p) the competition interface (p is 



enclosed by two semi-infinite maximizing polymers 
and starting from 7q and 7q, respectively, and with 
the same inclination Q. Therefore x < ^ in this case. 

Conclusions The connections studied above between 
the competition interface, the second class particle and 
maximal polymers fit into the interplay between the fiuc- 
tuation statistics and the global geometry of the growth 
interface developed by Prahofer and Spohn If the 

final macroscopic profile is curved then the competition 
interface follows a random direction (characteristic) in- 
tersecting the final surface at a point with non-zero cur- 
vature. In this case we have the KPZ scaling and the 
competition interface gets the transversal fluctuations, 
indicating the exponent x = 2/3. If the macroscopic 
profile is not curved we have two different situations. In 
the fiat case (stationary growth) the competition inter- 
face also follows the characteristics of the associated hy- 
drodynamic PDE and we still have the KPZ scaling. In 
the shock case the competition interface gets the longi- 
tudinal fluctuations which, in this case, are produced by 
the Gaussian fluctuations (x = 1/2) of the initial pro- 
file. On microscopic grounds one might suggest different 
rules for growth and competition. By universality we 
expect that from the knowledge of the curvature of the 
final macroscopic shape one can infer the asymptotics of 
the competition interface. This fits with the exponents 
founded by Derrida and Dickman Q in the Eden con- 
text since, in this case, the macroscopic profile is curved 
for angles 6 > 180° (we notice that in their simulations 
they have considered periodic initial conditions and so 
the longitudinal fluctuations in the shock direction are 
governed by the exponent 1/3 ^3)- 
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